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ABSTRACT 

Research  on  fundamental  mechanisms  of  erosion  is  the  incentive  for  develop¬ 
ment  and  application  of  the  numerical  methods  discussed  here.  The  methods 
iifck  particulate  statistical  dynamics  to  continuum  turbulent  reactive  flow 
and  solid  material  dynamic  (erosive)  surface  response.  The  modeling  incor¬ 
porates  calculated  thermal,  chemical,  and  mechanical  dynamic  processes  in 
solid,  liquid,  gas,  and  mixed  phr^es.  The  calculations  match  the  multi- 
component,  turbulent,  chemically-reactive  wall  boundary  layer  with  the 
unsteady  eroding  interface  boundary  using  microscopic  transport,  accommoda¬ 
tion  and  energy  distribution  sub-models.  Gas-borne  particle  motions  are 
coupled  to  the  turbulent  fljw  in  the  inertial  core  and  to  the  dissipative 
flow  near  the  wall  surface  using  two  basically  distinct  numerical  methods! 

(1)  particulate  turbulent  gas  flow,  statistical  dynamics  and  trajectory 
determination  in  the  inertial  core  region;  (2)  solution  of  the  detailed 
boundary  layer  heat  and  mass  interface  transport  with  boundary  conditions 
and  exchange  properties  modified  for  the  unsteady  influence  of  particle 
distributions  and  probabilistic  surface  dynamic  response.  Force,  energy 
and  motion  coupling  between  particles  and  gas  ie  assumed  weak,  with  negli¬ 
gible  back-influence  from  particles  to  gas  flow  structure.  The  present 
analysis  is  thereby  restricted  to  relatively  low  mass  loading  which  may  be 
found  in  the  ullage  region  between  propellant  bed  and  bounding  walls, 
propellant  face  and  projectile  base  as  well  as  between  propellant  base  and 
breech.  Particle-to-gas  relative  f luctuational  velocities  are  obtained  at 
incremental  time  steps  by  numerical  randomization.  They  are  subject  to  the 
constraint  that  their  ensemble  averaged  motions  are,  in  sum,  equivalent  to 
the  locally,  apriori  computed  mean  turbulent  intensity.  Results  also 
include  some  preliminary  molecular  solid  lattice  statistical  dynamics 
simulations.  These  address  detailed  estimation  of  the  microscopic  momentum 
and  energy  deposition  and  partition  between  a  rapidly  heated  surface' sub¬ 
jected  to  coincident  particle  collisions. 

INTRODUCTION  /(/ 

Considerable  attention  has  been  paid  to  recently  improved  numerical  Simula-  ^ 
tions  of  gun  barrel  interior  ballistics.  Among  these,  promising  detailed 
averaging  techniques  have  been  developed  for  computing  two— phase  flow. 
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etrese,  tersperature  and  bulk  thermochemistry  that  appear  to  dominate,  at 
le&at  ia  the  large,  the  one-dimensions!  propellant  bed  combustion— convec— 
tion  and  condensed  phase  motions,  interior  pressure  development,  and 
projectile  accelerations#  Experimental  correlation  of  some  aver- 

ege  predictions,  particularly  interior  pressures  and  projectile  motions, 
help  to  confirm  the  potential  of  these  methods  as  tools  for  gun  technology 
improvement.5  Use  of  these  promising  results  as  a  data  base  for  develop¬ 
ing  initial  and  boundary  conditions  for  subsequent  detailed  modeling  of 
erosive  flow,  transport  and  wall  material  response  appears  reasonable  and, 
to  sose  extent,  reassuring.  Average  dynamics,  equation  of  state  proper¬ 
ties,  burning  rates,  and  bulk  or  volume  averaged  stress  and  temperature 
predictions  are  available  for  checking  the  consistency  of  detailed  trans¬ 
port  and  kinetics  calculations  and  predictions. 

One  of  the  considerations  in  propellant  erosion  studies  is  the  influence  of 
the  populations  of  gas-borne  particulates  in  the  ullage  regions  surrounding 
the  combusting  propellant  bed  radially  and  the  regions  just  ahead  of  and 
just  behind  the  bed,  Fig.  1.  It  seems  intuitively  reasonable  that  unbumed 
(or  burning)  solid  propellant  particulates  impinging  on  or  abrasively 
sweeping  the  barrel  walls  will  enhance  wall  erosion  by  augmenting  heating 
and  shear  and  exothermic  gas-surface  chemical  processes  (by  continuous 
exposure  of  surface  reaction  sites).  It  has  been  observed,  hbwever,  that 
erosion  even  in  the  presence  of  these  "hostile"  particulates  can  be  coun¬ 
tered  to  some  extent  by  certain  "benign"  particulate  additives  tjhich  act  in 
some,  as  yet,  incompletely  understood  manner  (momentum  transports  buffer¬ 
ing?,  thermal  sinks?,  thermal  radiation  buffers?,  reduction  of  transport 
gradients?,  paticulate  cloud  chemical  absorption  layers?,  particulate 
Dufour  influences?,  turbulent  energy  and  vorticity  dissipation 
nuclei?). 


The  mechanisms  inherent  in  particulate  augmentation  of  erosion  and  additive 
particulate  erosion  suppression  are  the  central  issues  in  our  current 
theoretical  analysis.  We  abstract  the  problem  to  one  of  determining  solu¬ 
tions  of  the  coupled  particulate  laden  flow  and  transport  in  the  ullage 
region  between  propellant  bed  and  projectile.  We  further  simplify  initial 
and  boundary  conditions  at  a  planar  non-moving  porous  interface  with  flux, 
temperature  and  pressure  histories  imposed  so  as  to  be  consistent  with  that 
predicted  by  two-phase  NOVA  code  internal  ballistics  solutions.3*4 

The  gas  phase  turbulent  field  generated  by  such  a  propellant  flow  is 
modeled  as  that  which  would  be  initiated  by  a  porous  multiple  jet  plate 
turbulence  generator10  decaying  axially  as  implied  by  the  experiments  and 
by  a  conceptually  similar  turbulent  wake  decay  analysis.11  The  initial 
and  boundary  conditions  of  this  idealized  abstraction  are  summarized  in 
Fig,  2. 

PROCEDURE 

MODEL  FLOW  CONDITIONS 

In  previous  studies12*13  emphasis  was  placed  on  numerically  simulating 
the  complete,  multi-dimensional  Navier-Stokes  flow  field  in  the  wake  of  an 
accelerating  projectile.  However,  for  our  current  purposes,  it  is  appealing 
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U(t)  =  U(t)  +  U"(t") 
t"<t 


Fig.  1  Physical  geometry  of  flow  region  examined. 
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Fig.  2  Initial  conditions  for  model  problems. 

to  restrict  our  attention  to  details  of  the  particle  dispersal  associated 
with  an  experimentally  verified  pure  turbulent  flow.  Hence  we  remove 
consideration  of  the  other  complicating  flow  features  associated  with 
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accelerating  projectile,  projectile  wake,  and  propellant  bed  combustion 
irregular it ics . 14 

The  fluid  dynamic  situation  which  develops  includes  the  solid  aerosol 
consisting  of  sparse  distributions  of  urburned  propellants  and  erosion 
•uppreesive  additives  such  as  finely  divided  Ti02,  Si02,  and  talcum 
powder. 


Particulate  loading  considered  in  the  present  analysis  is  that  associated 
with  sparse  population  densities  of  the  order  of  103  cm-3  or  less. 

These  relatively  low  mass  loading  levels  together  with  characteristic  aero¬ 
sol  sire  selected  (0.1  ym  to  10  ym)  permit  treatment  of  weak  coupling  from 
the  carrier  gas  to  the  particle  motions  with  neglect  of  back  influence  of 
the  particles  on  the  carrier  gas.  The  original  gun  barrel  situation  is 
illustrated  in  Fig.  1  with  decomposition  of  the  flow  into  mean  (ft)  and 
fluctuating  (U"  *  u*)  components  as  noted.  Our  attention  is  on  a  simple 
conceptual  region  for  the  forward  ullage  region  (shown  cross-hatched)  at 
early  times.  Boundary  and  initial  conditions,  shown  in  Fig.  2  include  an 
accelerating  inertial  mean  axial  flow  (Us),  a  gradual  cooling  of  the 
accelerating  propellant  products  as  implied  by  the  time  derivative  of  the 
interior  temperature  (T)  from  an  initial  flame  temperature  of  3200  K,  the 
ia  initially  cool  (T  ■  300  K) ,  entrance  conditions  are  for  a  plane 
adiabatic  duct  (Q  “  heat  flux)  and  the  turbulence  as  represented  by  the 
kinetic  energy  (k)  at  the  jet  plate  origin  is  constant  and  homogeneous 
(immediately  well-mixed). 

COKTIKUUM 

The  basic  Navier-Stokes  solutions  subject  to  initial  and  boundary  condi¬ 
tions  previously  described  are  obtained  using:  (1)  a  combination  of  the 
explicit  MacCormack  differencing  scheme  together  with  the  two  equation 
model  closure  of  Wilcox  and  Traci  and  (2)  the  Briley-McDonald  factored 
implicit  scheme  with  a  Jones— Launder  two  equation  closure  model.  These 
produce  the  Navier-Stokes  field  quantities,  together  with  turbulence 
kinetic  energy,  ic,  and  dissipation  rate,  £.12,13 

Mass  conservation  considerations  for  turbulent  flow  require  global  continu- 
ity,  eq.  (1),  as  well  as  both  scalar  time  averaging,  eq.  (2)  and  vector 
velocity  mass  averaging,  eq.  (3).  In  the  boundary  layer  region,  gas  phase, 
gas-surface,  and  surface  chemical  reactions  require  additional  species 
conservation  and  production  (ii>e  f  0)  as  well  as  multicomponent  species 
diffusive  exchange  considerations,  eq.  (4). 

Global 
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Substitution  of  the  velocity  perturbations  in  the  compressible  Navier- 
Stokes  equations,  scalar  averaging  (eq.  (2))  and  mass  averaging  (eq.  (3)) 
and  applying  global  continuity  (eq.  (1))  results  in  the  compressible 
Reynolds  equation  for  time-averaged  momentum, 

Momentum 

a(pUi>+  d 

at  +  dXj 


In  addition  to  the  Reynolds’  stress  terra,  (Xj’-jUj  *  R£ j ,  eq.  (5) 
introduces  void  fraction  weighting,  defined  byeq.  (6),  and  developed  by 
writing,  separately,  the  equations  for  particle  motions  and  gas  phase 
motions,  then  combining  them  through  simplified  force  and  condensed  phase 
production  (weak  coupling)  considerations  to  form  the  mixture  representa¬ 
tion,  eq.  (5)..  To  complete  the  weak  coupling  between  turbulent  field  and 
the  particle  field  requires  that  an  explicit  evaluation  must  be  made  of  the 
forces  on  the  particles  induced  by  the  turbulent  motions.  This  will  be 
addressed  subsequently.  At  this  stage  one  may  note  that  the  stress  terms 
in  the  elementary  mixture  balance  equation  are  simply  weighted  by  the  void 
fraction.  The  weighting  gives  the  correct  vapor  phase  limit  as  the  rela¬ 
tive  volume  occupied  by  the  particles  becomes  negligible  (Vp/Vg%  0) 
and  implies,  correctly,  as  discussed  by  Crowe14,  that  for  light  particle 
mass  loading  the  particles  contribute  negligibly  to  the  mixture  stress 
tensor.  The  coupling  is  weak.  There  is,  as  anticipated,  no  back  influence 
develops  from  the  particles  to  the  mixture. 

In  terms  of  the  tctal  enthalpy,  H  (the  sum  of  internal  energy,  pressure 
work,  and  mean  flow  kinetic  energy)  the  corresponding  energy  equation  for 
the  turbulent  mixture  of  particles  and  vapor  is, 


i  Uj  Uj  +  ©Sj,  p — 0  (fjj  -  p  ’ 
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Figure  3  illustrates  the  results  of  the  continuum  turbulent  Navier  Stokes 
solutions  showing  qualitative  comparison  between  the  calculated  results  and 
Klingenberg's15  experiment  in  bore  measurements.  Figure  4  illustrates  a 
similar  quantitative  comparison  for  the  temperature  field.  Both  flow  cal¬ 
culations  are  taken  from  earlier  results  which  include  complicating  fea¬ 
tures  such  as  the  projectile  turbulent  wake.12*13  These  vortical  inter¬ 
action  features  have  been  purposely  omitted  from  our  current  analyses  which 
use  the  simplified  grid  plate  induced  wind  tunnel  turbulence  to  isolate 
particulate  turbulent  coupling  and  transport.16  Radial  distributions  of 
the  time  dependent  turbulent  intensity  are  illustrated  in  Fie.  5  and 
compared  to  the  quasi-steady  jet  mixing  experiments  of  Owen.*7  The 
Jones-Launder  (J-L)  model  closure  relations  appeared  to  give  qualitatively 
closer  agreement  to  the  experiments  than  did  the  Saf fman-Wilcox-Traci 
(S-W-T)  results.  The  Owen  results  are  not  appropriate  for  comparison  with 
the  near  centerline  predictions,  since  the  experimental  set-up  used  an 
isolated  round  jet  as  opposed  to  the  homogeneous  multiple 


(G.  K I  ingen  berg  and  H.  Mach,  1978) 

Fig.  3  Comparison  of  predicted  and  experimental  velocity  profiles. 
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Fig.  4  Comparison  of  predicted  and  experimental  temperature  profiles. 


Fig.  5  Radial  turbulent  intensity  profiles. 

jet  mixing  region  found  near  the  turbulent  grid  plate.  For  comparison,  the 
core  and  near  centerline  distributions  the  Betchov-Lorenzen^  centerline 
axial  distributions  from  their  experiments  were  plotted  against  our  theor¬ 
etical  predictions.  Good  agreement  is  found  at  several  axial  positions 
(measured  from  the  grid  plate  origin  in  the  downstream  direction).  These 
comparisons  were  reported  in  an  earlier  summary. 
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BOUNDARY  LAYER  MATCHING 


The  wall  boundary  region  is  dominated  by  the  steepest  fluid-thermodynamic 
gradients  yet  its  extent  is  small  in  comparison  to  the  overall  bore  dimen¬ 
sions.  In  the  current  approach  the  multi-component  reactive  turbulent 
boundary  layer  solutions  with  wall  ablation,  recession,  blowing  and  melt 
erosion  are  obtained  by  applying  a  separate  non-similar,  reactive  boundary 
layer  method.  Boundary  conditions  are  supplied  by  the  core  flow  continuum 
solutions,  together  with  the  computed  statistical  particle  concentration 
and  energy  distributions.12  The  matching  of  the  separate  computational 
domains,  each  with  its  own  distinct  formulation  and  me'. hod  counters  the 
numerical  difficulties  induced  by  the  disparate  time  and  length  scales 
associated  with  each  region.  The  solutions  for  Navier-Stokes  reacting  core 
flow  and  boundary  layer  flow  are  asymptotically  matched  in  their  overlap 
domain.  Figure  6  illustrates  the  matching  procedure.  For  "level"  2  addi¬ 
tional  conditions  are  needed  and  found.  These  require  minimization  of 
and  together  with  the  wake  flow  velocity  and  thermodynamic  edge  condi¬ 
tions  to  locate  the  wall  boundary  layer/core  interface.12*1^  The  match¬ 
ing  procedure  depends  on  the  singular  perturbation  nature  of  the  core 
generated  turbulent  region  as  the  "outer  solution"  (the  circled  "e"  region) 
in  relation  to  the  wall  generated  turbulent  region  as  the  "inner  solution" 
(the  circled  'V  region),  the  overlap  region  to  which  the  asymptotic 
matching  is  applied  is  shown  cross-hatched. 

The  other  matching  "levels’*  shown  in  Fig.  6  include  level  (1)  the  ordinary 
(first  order)  boundary  layer  to  core  potential  flow  matching  and  level  (3) 
which  invokes  an  explicitly  computed  large  scale  turbulence  structure, 

Ue,  as  well  as  the  small  scale,  Ue  ,  modeled  as  a  continuum. 

Level  (3)  represents  a  subsequent  step  in  the  modeling  which  has  not  yet 
been  completed. 
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Fig.  6  Schematic  of  boundary  layer-core  overlap  procedure. 
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•  vS?ATIgT£CAL  1U1BULEHT  PARTICLE  MOTIONS 

'  &  £ea tree  of  the  present  numerical  simulations  is  the  coupling  of  an 

-££T&y  of  particles  to  the  Reynolds  averaged  Nsvier-Stokea  gas  phase 

The  statistical  probability  distributions  associated  with  predic¬ 
ted  localised  (space  and  time)  turbulent  kinetic  energy  and  dissipation 
rcte  results  are  obtained  by ‘ numerical  random  (Monte  Carlo)  sampling  of  the 
turbulent  excursion  velocities.^  The  random  perturbstional  velocity 
distribution  is  subject  to  a  constraint  that  the  ensemble  average  must 
correspond  to  the  separately  generated  continuum  flow  field  prediction  of 
local  turbulent  intensity* 

Aerosol  eetion©  are  treated  incrementally  in  time,  ^t,  following  selection 
of  the  excursion  velocities  from  the  distribution  sampling.  Particle 
forces  are  then  computed  dependent  on  these  statistical  excursion  veloci¬ 
ties.  The  particle  forces  derived  by  Crove*^,  are  rewritten  in  a 
Lsgrangian  framework  ©ovin;  with  the  mean  gas  dynamic  flow,  U£.  Forces 
considered  include  drag,  virtual  mass,  Basnt,  Magnus,  and  Saffman  Lift,  in 
order,  as  indicated  by  eq.  (13).  These  forces  drive  the  particles  in 
linear  motions  over  trajectories  limited,  automatically,  by  the  smallest 
time  increment  between  collisions.  Aggregate  incremental  motions  simulate 
particle  dispersion,  momenta,  energy,  and  wall  collisions.  Hard  sphere 
particle-particle  collisions  are  prescribed  for  particle  encounters  in  the 
computations. 


Wd3v^/’*l  (u,  _  Vjl  VV  + 


ffd2 

*«k  —  Pf ««;  -  °p.j>  wk  + 


3u? 


c2VupTd2{u;-«pJ} 


(13) 


The  random  sampling  is  applied  for  1100  particles  with  expectancies  distri¬ 
buted  among  110  bins.  The  random  velocity  at  any  point  ro  is  a  sum  of 
the  sampled  velocities,  <6u*>  plus  the  auto— correlation  to  the  original 
(tD)  velocity  fluctuation  for  the  "memory”  of  the  turbulent  eddy  as  it 
rvolves  over  the  sampling  time  interval,  At.  This  is  shown  in  eq.  (14) 
together  with  a  schematic  representing  the  bin  subdivision  over  the  samp¬ 
ling  interval  (A  =  At).  The  base  distribution  is  considered  to  be  the 


to 


normal  distribution,  eq.  (15).  Selective  bin  weighting  is  applied  using 
the  "importance  sampling"  techniques  developed  in  chemical  physics 
applications  to  selectively  sample  the  most  interactive  portions  of  the 
spectra,  hence  speed  convergence  without  exhaustively  large  numbers  of 
selections. Trial  weighting  functions  are  listed  in  eq.  (16).  Best 
results  for  the  present  cases  were  achieved  using  the  parabolic  bin 
weighting,  vj  *  2(1  - which  favors  the  population  near  the  origin 
(t0>‘ 


u*  (r0,t0  +  At)  »  u*  (r0,tQ)  R*uu  (rQ#  At)  +  <5u#> , 

fl(u-u)  ■  (Y6u#  yflvf'  x 

exp  (— 6u#2/2Y|u*)  , 

Bin  weight*: 


(14) 


(15) 


w, 


f  1. 

Zt/2, 
1/2-6-  > 

3(1  — fr), 
2(1 --fr2). 


(16) 


Specific  emphasis  of  these  calculations  is  on  description  of  particle 
concentration  size  and  mass  distributions  together  with  their  resulting 
infuences  on  local  flow  composition,  state  aid  transport  properties,  as 
well  as  on  gas-to-wall  surface  activity,  energy  and  momentum  transport. 

Figure  7  shows  the  variance  of  the  lateral  dispersion  about  assumed 
centerline  injection,  predicted  by  the  numerical  simulations  and  applied  to 
0.3  ym  and  5  ym  aerosol  particulates.  The  smaller  particles  relax  with 
relative  rapidity  to  the  local  turbulent  flow  motion  so  that  their  disper- 
sal  pattern  almost  matches  that  predicted  by  assuming  an  instantaneous 
equilibration  to  the  local  gas  velocities.  The  larger  particles,  however, 
lag  the  local  flow  currents,  creating  narrower  dispersion  patterns  and 
shortened  penetration  distances. 

Figure  8  illustrates  the  effects  of  the  initial  particle  distributions. 

From  among  the  example  packing  distributions  previously  analyzed,16  an 
initial  packing  conceptually  arranged  outside  -the  flaps  is  selected  as 
representative  of  one  cf  the  most  effective  configurations  tested  to 
achieve  maximum  dispersal  of  "benign"  (erosion-suppression)  additives. 

Size  effects  are  evident  in  assessing  the  accumulated  percentage  of  the 
original  particle  population  reaching  the  wall. 
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Fig.  8  Turbulent  flow  particle  dispersal  for  indicated  initial 
distribution. 


BOUNDARY  LAYER  TRANSPORT  AND  WALL  ACCOMMODATION 


Additive  particulate  or  "hostile”  particulate  concentrations  are  incorpor¬ 
ated  in  the  edge  conditions  through  changes  in  the  edge  entropy  profile, 
the  corresponding  gas  state,  and  the  altered  conditions  on  the  momentum 
integral  due  to  increased  mass  flux.  Within  the  boundary  layer  the  trans¬ 
port  exchange  coefficients  for  multi-component  species  are  adjusted  with 
"frozen"  contributions  from  the  non-equilibrium  thermal  and  momentum  relax¬ 
ation  of  the  particles.  Figure  9  shows  the  effects  of  two  sizes  of  addi¬ 
tive  particles  on  boundary  layer  enthalpy  profiles  and  on  the  wall  shear. 
The  indications  are  that  twenty-five  fold  decrease  in  particle  size 
corresponds  to  a  15%  to  20%  decrease  in  erosive  wall  heat,  shear  and  mass 
transfer.  In  the  absence  of  any  particles  the  corresponding  levels  in 
erosive  transport  increase  up  to  40%. 

A  more  detailed  exposition  of  the  particle  influences  on  the  transport 
processes,  together  with  the  theoretical  perturbation  method  used  to 
develop  the  particle-augmented  gas  transport  properties  is  currently  in 
preparation  for  a  subsequent  paper. ^  In  addition,  further  illustrations 
of  the  effects  of  coincident  wall  surface  micro-response  and  erosion  on  the 
transport  processes  will  be  provided  in  this  subsequent  paper.  Here  we 
will  limit  our  discussion  to  some  preliminary  description  of  the  modeling 
of  additive  laden  gas  to  surface  accommodation  and  energy  partition. 

The  wall  surface  accommodation  of  energy  and  momentum  associated  with  the 
arrival  of  the  particulate  cloud  (a  continuum  distribution  of  particulates) 
can  be  modeled  using  a  number  of  more  or  less  successful  engineering  models 
on  gas  to  surface  transport.  In  the  present  analysis  however,  adoption  of 
a  particular  modeled  process,  a  priori,  is  self-defeating  since  the  funda¬ 
mental  erosive  mechanisms  under  investigation  are  imbedded  in  the  model  and 
its  governing  assumptions. 

To  circumvent  over ly-restr icuive  assumptions  and  to  help  unfold  the  funda¬ 
mentals  of  gas-particulate-surface  interaction,  current  use  is  made  of 
molecular  dynamics  lattice  collision  calculations.  These  first  principle, 
ab  initio,  molecular  dynamics  calculations  are  well-established  procedures 
developed  and  refined  by  Karo,  et  al.*^  at  Lawrence  Livermore  Laboratory 
to  investigate  details  of  kinetics,  detonation,  and  combustion  on  a 
molecular  level. 

To  supplement  our  knowledge  of  particulate  laden  gas  surface  energetics, 
Karo  is  applying  this  numerical  statistical  dynamic  technique  to  the 
description  of  particulate  collisions  and  energy  deposition  with  modeled 
(brittle)  wall  surface  layers  of  atoms  simultaneously  heated  by  boundary 
layer  transport.  To  represent  the  thermally  excited  surface  the  lattice  (a 
community  of  particles  connected  by  modeled  intermolecular  potentials)  is 
set  into  vibrational  motions.  Resonances  induced  by  particulate  collisions 
(the  "particle"  is  another  colliding  lattice)  are  followed  incrementally  in 
time  during  the  calculation.  Energy  deposition  distributions,  stress 
development  and  dynamics  of  lattice  deformation  which  may  signity  failure 
are  revealed  as  point  defects  (spall,  dislocations,  etc.)  or  as  averages 
(in  the  ensemble  sense)  which  may  be  extended  to  a  macroscopic  model  cf 
plastic  deformation,  failure  and  surface  removal. 
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Figure  10  is  an  illustration  of  an  impact  simulation  showing  the  collision 
of  a  particle  at  the  left  hand  side  of  a  lattice  group  representing  a 
surface  layer  of  8  molecules.  The  surface  particle  response  to  collision 
motions  are  frozen  in  the  figure  at  a  lattice-related  characteristic 
frequency  period  of  39.  The  information  generated  by  a  systematic  study  of 
the  time  histories  simulating  a  variety  of  such  collisions  and  surface 
layer  models  are  used  to  develop  energy  accommodation  and  partition 
information  for  the  further  development  and  modeling  of  erosive  inducing 
transport  mechanisms  at  the  solid  wall  surface. 

SUMMARY  AND  CONCLUSIONS 

The  present  discussion  emphasizes  the  matching  of  statistical  mechanical 
and  continuum  modeled  processes  to  help  assess  the  role  of  particulates 
(both  hostile  and  benign)  on  gun  surface  wall  erosion.  Particulate 
influences  on  erosion  is  offered  here.  The  calculation  and  modeling  at 
this  time  have  emphasized: 

o  Effects  of  turbulence  on  species  concentrations)  thermal  spatial 
distributions,  and  non-equilibrium  relaxation  of  the  particle 
field; 

o  Transport  exchange  influences  of  particulate  clouds  in  the 
turbulent  eroding  wall  boundary  layer  region; 

o  Modifications  to  properties  for  simultaneous  heat,  mass  and 
particulate  transfer,  energy  accommodation  and  partition  on 
brittle  surface  layers  based,  in  part,  on  gas-solid  statistical 
dynamics  considerations. 

The  primary  (erosive-suppressive)  influence  of  benign  (non-burning) 
particles,  if  small  enough  to  disperse  rapidly  with  significant  populations 
to  the  walls  appears  to  be  in: 

o  Mass-inertial  and  volumetric  constraints  act  to  reduce  the 

molecular  transport  processes  associated  with  the  gas  field; 

o  Particle  flow  in  the  boundary  layer,  comprising  a  moving  curtain 
scatters  incident  particulates  and  fragments,  preventing  their 
reaching  the  wall; 

o  Particle  curtain,  associated  with  small  diameter  particles, 

offers  considerably  greater  surface  area  than  does  the  exposed 
wall;  hence  is  an  effective  chemical  absorption  layer  for  hostile 
exothermic  reactions,  in  addition  it  provides  a  scattering 
barrier  for  incident  thermal  radiation,  effectively  shielding  the 
wall  from  these  hostile  influences; 

o  Particles  in  the  core  flow  provide  acoustical  damping  for  the 

high  frequency  pressure  oscillations  which  drive  the  fluctuating 
velocity  field,  hence  act  to  damp  (lower  the  intensity)  of  the 
combustion  originated  turbulence,  and  thereby  lower  the  eddy 
transport  of  erosive  heat  and  shear  effects  on  wall  erosion. 
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NOMENCLATURE 

Modified  Stokes  drag  coefficient  on  particles 
Effective  particle  sphere  diameter 
ith  component  of  particle  force 

Probability  density,  frequency  function  in  statistical 
weighting 

Integral  average  of  a  random  property  distribution 

Particle  internal  energy 

Number  of  samples  in  statistical  average 

Temperature 
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rj 


*1 
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ic 
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Jth  cespoaent  of  mean  velocity 
jfch  costpeneat  velocity  (fluctuation) 

floaty,  «*  floctmtin* 

jth  cof^oaent  vector  distance 
Eeysolda  stress  tensor 
ith  random  variable  in  •  turbulent  field 
Turbulence  dissipation  rate 
Arbitrary  fluctuating  component 
Turbulence  kinetic  energy 

Void  fraction,  ratio  of  gas  volume  to  mature  volume 
Geo  molecular  viscosity  coefficient 
Cfio  density 

Particle  rotational  velocity  component 
VclU222 


